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26 


Abstract 


27 In this study, MODIS fine mode fraction and MISR non-spherical fraction are 

28 used to derive dust and smoke AOT components (idust and x sm oke) over the tropical 

29 Atlantic, and their variabilies related to the Madden- Julian Oscillation (MJO) are then 

30 investigated. 

3 1 Both MODIS and MISR show a very similar dust and smoke winter climatology. 

32 Xdust is found to be the dominant aerosol component over the tropical Atlantic while x smo ke 

33 is significantly smaller than Xdust- The daily MODIS and MISR Xdust are overall highly 

34 correlated, with the correlation coefficients typically about 0.7 over the North Atlantic. 

35 The consistency between the MODIS and MISR dust and smoke aerosol climatology and 

36 daily variations give us confidence to use these two data sets to investigate their relative 

37 contributions to the total AOT variation associated with the MJO. However, unlike the 

38 MISR dust discrimination, which is based on particle shape retrievals, the smoke 

39 discrimination is less certain, based on assumed partitioning of maritime aerosol for both 

40 MISR and MODIS. 

41 The temporal evolution and spatial patterns of the Xdust anomalies associated with 

42 the MJO are consistent between MODIS and MISR. The Xd us t anomalies are very similar 

43 to those of x anomalies, and are of comparable magnitude. In contrast, the M JO-re la ted 

44 x smo ke anomalies are rather small, and the x mar anomalies are negligible. The consistency 

45 between the MODIS and MISR results suggests that dust aerosol is the dominant 

46 component on the intra-seasonal time scale over the tropical Atlantic Ocean. 
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47 1 . Introduction 


48 The Madden- Julian Oscillation (MJO) [Madden and Julian, 1971; 1972] is the 

49 dominant form of the intra-seasonal (30-90 day) variability in the tropical atmosphere. It 

50 is characterized by slow (~5 m s' 1 ) eastward-propagating, large-scale oscillations in the 

51 tropical deep convection over the equatorial Indian Ocean and western Pacific during 

52 boreal winter (November-April) [Lau and Waliser, 2005; Zhang, 2005]. Recently, there 

53 is an emerging strong interest in the impacts of the MJO on atmospheric composition 

54 [Tian and Waliser, 2011], such as aerosol [Tian et al., 2008; 2011], ozone [Tian et al., 

55 2007; Weare, 2010; Li et al., 2011], carbon dioxide [Li et al., 2010], and carbon 

56 monoxide [Wong and Dessler, 2007]. 

57 Tian et al. [2008] first examined the aerosol variability related to the MJO using 

58 global aerosol products from multiple sensors on various satellite platforms. That study 

59 revealed large intra-seasonal variations in the satellite-derived aerosol products over the 

60 Tropics, though the interpretation in tenns of actual aerosol behavior was ambiguous. 

61 Tian et al. [2011] further investigated the MJO-related aerosol variability over the 

62 tropical Atlantic Ocean using the aerosol optical thickness (AOT) product from the 

63 MODIS (Moderate Resolution Imaging Spectroradiometer) on the Aqua satellite. They 

64 suggested that the MJO-related intra-seasonal variance accounts for about 25% of the 

65 total AOT variance over the tropical Atlantic. They also found that the AOT anomalies 

66 are negatively correlated with the low-level zonal wind anomalies over most parts of 

67 tropical Atlantic, with the latter leading the fonner by about 6 days. This indicates the 

68 MJO may modulate Atlantic aerosol transport through its influence on the Atlantic low- 

69 level zonal winds. Given the potential predictability of the MJO extending to 2-4 weeks 
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70 [e.g., Waliser, 2005], the study by Turn et al. [2011] implies that the Atlantic aerosol 

7 1 concentration may be predictable with lead times of 2-4 weeks, which in turn may lend 

72 important guidance to predicting air quality, dust storm activity, and ocean nutrient 

73 deposition over the Atlantic Ocean. 

74 Nevertheless, Tian et al. [2011] examined only the total AOT anomalies and did 

75 not consider the contribution of different aerosol types to the total AOT anomalies. It is 

76 well known that the aerosol over the tropical Atlantic Ocean in the boreal winter season 

77 is a mixture of mineral dust from the Sahara desert and the Sahel region, biomass burning 

78 smoke from the Sahel and African savanna regions, and marine aerosol (primarily sea salt 

79 and secondly sulfate aerosols) from the ocean surface [Kaufman et al, 2002; 2005a; 

80 2005b]. Since dust, biomass burning smoke and marine aerosols play very different roles 

81 in the radiative forcing and cloud fonnation process, it is of great interest to partition the 

82 total aerosol into individual aerosol components and examine the MJO-related variability 

83 in each aerosol type. 

84 Previous studies [Kaufman et ah, 2002; 2005a; 2005b] have suggested that 

85 satellite data distinguishing fine -mode aerosols from coarse-mode aerosols could be used 

86 to separate the aerosol into specific types, since different aerosol types (e.g., smoke, dust, 

87 and maritime sea salt aerosols) have different fine mode fraction (FMF) values. FMF is 

88 the fraction of total AOT contributed by the fine-mode aerosols. For example, smoke 

89 from wildfire biomass burning and urban-industrial air pollution from human combustion 

90 processes are dominated by fine-mode particles. On the other hand, desert dust and sea 

91 salt, mainly a result of wind and erosion processes, have a much larger proportion of 

92 coarse-mode particles. Kaufman et al. [2002; 2005a; 2005b] first demonstrated the 
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93 usefulness of combining the total AOT and FMF measurements from MODIS/Terra to 

94 estimate dust and biomass burning smoke AOT components (see detailed description in 

95 section 2). This method has since been widely adopted by the community to understand 

96 aerosol types and their climate forcing [ Bellouin et al., 2005, among many others]. 

97 Measuring the reflectance of a target from different directions is very useful 

98 because geophysical media, for example aerosols, reflect solar light differently in 

99 different directions. The variations between the reflectances acquired from a variety of 

100 observation angles can be interpreted (with appropriate models) in tenns of aerosol 

101 properties such as particle size, shape, and single-scattering albedo [Kahn et al., 1998, 

102 2001; Chen et al. , 2008]. In particular, MISR’s sensitivity to the characteristics of the 

103 aerosol scattering phase functions enables it to distinguish between the non-spherical and 

104 spherical particles, and thus provides a possible way to separate mineral dust aerosols 

105 from other aerosol components. A series of studies has explored the ability of MISR to 

106 retrieve mineral dust properties theoretically [Kahn et al., 1997, 2001; Kalashnikova et 

107 al., 2005; Kalashnikova and Kahn, 2006], the sensitivity of the theoretical results, as well 

108 as the application of non-spherical dust models for five Saharan dust field events over the 

109 Atlantic Ocean [ Kalashnikova and Kahn, 2008], 

110 In this study, we use MODIS FMF to separate the total column AOT into its 

111 individual types over the tropical Atlantic, especially the dust and biomass burning 

112 smoke aerosols, because of its more frequent sampling. We will also use MISR aerosol 

113 non-spherical fraction to derive the dust aerosol component because of its greater particle 

114 type information content, and compare it with that derived from MODIS. The resulting 

115 intra-seasonal variabilities associated with the MJO in dust and biomass burning smoke 
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116 aerosols over the tropical Atlantic Ocean are then examined. 

117 The rest of this paper is organized as follows. Section 2 describes the MODIS and 

118 MISR data, and the methodology used to derive specific aerosol components from the 

119 total AOT. The climatology of specific aerosol component as well as the comparison of 

120 dust aerosol between MODIS and MISR is presented in Section 3 in order to examine the 

121 fidelity of the methods described in Section 2. The main results of this paper, the MJO- 

122 related dust and smoke aerosol anomalies, are presented in Section 4. Conclusions and 

123 discussions are presented in Section 5. 

124 2. Methodology and Data Description 

125 2.1. MODIS 

126 2.1.1. Method to Calculate Dust, Smoke AOT 

127 The method developed by Kaufman et al. (2005a, b) to derive the AOTs for dust 

128 and biomass burning smoke using MODIS total AOT and FMF measurements is briefly 

129 described as below. Over the tropical Atlantic Ocean during the boreal winter season, the 

130 aerosol is a mixture of dust, biomass burning smoke, and marine aerosols. Note that the 

131 term ‘anthropogenic aerosol’ in Kaufman et aids papers is replaced by ‘biomass burning 

132 smoke’ as the latter is far more appropriate tenn for the region of study. However, the 

133 ‘smoke’ does include a small contribution from air pollution that originates primarily 

134 from the US and South American continents. With the two constraints that both the total 

135 AOT and its fine model fraction can be partitioned into contributions from the dust, 

136 smoke, and marine aerosol components, we have the following equations: 

137 X T dust + Tsmoke "I" T mar (1) 

138 f x T — fdust X- Oust fsmoke X Umoke fmar X ^mar (2), 
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139 in which x and f denote AOT and FMF respectively (x and AOT as well as f and FMF are 

140 interchangeable in this paper), and the subscripts ‘dust’, ‘smoke’, and ‘mar’ indicate dust, 

141 smoke, and marine aerosol component, respectively. Rewriting (1) and (2), we get 

142 Xdust — [^^^(fsmoke - f) — T'mar x (f'smoke — fmar)]/(fsmoke — fdust) (3) 

143 Xsmoke — [X x (f—fdust) —, £mar X (fmar — fdust)]/(fsmoke — fdust) (4). 

144 With x and f being MODIS measurements, Xdust and x smo ke can be computed directly if fd ust , 

145 fmokc? friar, and x mar are known. 

146 In Kaufman et aids study, the FMFs of dust, smoke, and marine aerosols are 

147 assumed to be constant, and were derived by averaging the MODIS/Terra FMF 

148 measurements over selected regions and time periods where one specific aerosol type 

149 dominates, with their uncertainties estimated from these selected measurements: 

150 fd ust =0.5±0.05, f smo k e =0.9±0.05, and f mar = 0.3±0.1 [2005b]. However, the actual FMF 

151 values vary with season and location [Mating et al., 2003; Jones and Christopher, 2007, 

152 2011; Yu et al., 2009], thus large uncertainties are expected from using constant f dust , 

153 fsmoke, fmar when applying Kaufman’s formula. More discussion on this issue as well as 

154 the sensitivity of our results to the FMF values will be given in Section 4.3. 

155 The marine AOT, x mar , depends strongly on surface wind speed as its primary 

156 component is sea-spray salt [e.g., Smirnov et al., 2003]; it is estimated using the empirical 

157 formula in Kaufman et al. [2005b]: 

158 x mar = 0.007W+0.02 (5) 

159 Here, W is the surface (10 meter) wind speed from the ECMWF ERA-Interim reanalysis 

160 [ Dee et al., 2011], The global mean of x mar is around 0.06±0.005 [e.g., Kaufman et al.. 
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161 2001], 

162 2.1.2. MODIS Data Description 

163 MODIS instruments are operating on both the Terra and Aqua satellites. Each has 

164 a viewing swath width of 2,330 km and views the entire surface of the Earth every one to 

165 two days. In this study, daily AOT and FMF measurements at 0.55 pm from the 

166 MODIS/Aqua Level-3 Collection 5.1 dataset [Levy et al., 2009] on 1° x 1° spatial grids 

167 are used. The uncertainties of x are ± (0.03 + 0.05x) over ocean and ± (0.05 + 0.1 5x) over 

168 land. In MODIS, the fine-mode aerosols refer to the aerosols with a size distribution of 

169 radii centered between 0.1 and 0.25 pm, whereas the coarse-mode aerosols refer to those 

170 with a size distribution of radii centered between 1 and 2.5 pm. Validation studies 

171 indicate that, over ocean, the uncertainties of f are large for low AOT (x<0.15) but 

172 typically less than about 20% for large AOT (x>0.15) [ Kleidman et al., 2005; Renter et 

173 al., 2005], Over land, MODIS does not provide any quantitative information about the 

174 aerosol size [Levy et al., 2010]. As a result, we use the x and f data over the ocean only. 

175 The period of 4 July 2002 to 1 June 2009 is used for consistency with the study by Tian 

176 et al. [2011], 

111 Note that we use MODIS/Aqua aerosol data rather than MODIS/Terra data even 

178 though MISR is on board Terra. One reason is that we want to be consistent with the 

179 study in Tian et al. [2011], More importantly, MODIS/Terra and MISR aerosol retrievals 

180 are much less spatially overlapped due to the exclusion of sun glint regions over the 

181 ocean in MODIS. As a result there are many fewer days on which both MODIS/Terra and 

182 MISR made observations of a given region, compared to the situation between 

183 MODIS/Aqua and MISR. Since the direct comparison between daily MODIS and MISR 
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184 aerosol data serves as an essential part of this study, we chose to use MODIS/Aqua data. 

185 The results based on MODIS/Terra were compared with those based on MODIS/Aqua, 

186 and it is found that the difference is negligible for our application. 

187 2.1.3. MODIS Data Rejection 

188 In this work, some MODIS x and f measurements are rejected in our calculations 

189 for the reasons described below. 

190 The dependence of idust and x smo ke on x and f based on equations (3) and (4) is 

191 illustrated in Figure 1, where x mar is set to 0.06 (the approximate global mean value of 

192 x mar ). It is seen that the larger the x, the larger Xdust and x smo ke. Also, as f increases, x smo k e 

193 increases and Xdust decreases. Figure 1 also suggests that equations (3) and (4) can produce 

194 reasonable (nonnegative) values for both Xdust and x smo ke only when the paired x and f 

195 measurements fall within a limited region bounded by the white dashed lines (xd us t =_ 0.03 

196 and x smo ke =_ 0.03). We relax the limit of Xdust or x smo ke to -0.03 instead of 0 because the 

197 uncertainties of the MODIS AOT are ±0.03 over ocean, thus a small negative AOT up to 

198 -0.03 is regarded as indistinguishable from the value 0 [Remer et al., 2005], Note that 

199 calculated Xdust and x smo k e in the range of 0 to -0.03 have been set to zero. Outside this 

200 region, either Xd us t or x smo k e is too negative, and correspondingly, x smo k e or Xd us t would be 

201 larger than the total AOT, which is not physical. Therefore, the x and f measurements 

202 giving rise to such Xdust and x smo k e values are rejected. In addition, x measurements greater 

203 than 2 are excluded, given the large uncertainties from cloud contamination [Zhang et al., 

204 2005], 

205 The count distribution of paired MODIS x and f measurements as a function of x 

206 and f is summarized for the tropical Atlantic (20°S - 30°N, 60°W - 20°E) and for the 
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207 2002-2009 boreal winters (October to April) in Figure 2a. It is found that the majority of 

208 observations fall within the region bounded by two lines (idust = “0.03 and x smo ke = -0.03); 

209 however, there is also a large number of observations falling outside the two lines, which 

210 are rejected in the calculation. The rejection results from both the limitations in MODIS 

211 data and the limitations in Kaufman’s method. For example, in some cases x is typically 

212 between 0. 1-0.3, which is much larger than the baseline x mar , while f is extremely small, 

213 nearly 0 (see bottom of Figure 2a). These observations are very likely artifacts of cloud 

214 contamination [e.g., Zhang et al., 2005; Tian et a/., 2008], For other cases, x is also 

2 1 5 significantly larger than x mar and f falls between f mar and fd us t, suggesting the aerosol is a 

216 mixture of dust and sea salt (see the area between the line of x smo ke = “0.03 and the line of 

217 f mar ). These are likely valid observations; however, due to the difficulties in separating 

218 dust from sea salt over this regime [ Kaufman et al. 2005b], we don’t include these 

219 observations in our calculations. 

220 The percentage of days with MODIS x and f measurements available relative to 

221 the total days during 2002-2009 boreal winters (1269) is shown in Figure 2b. The 

222 occurrence of cloudy conditions affecting the retrieval of MODIS aerosol parameters is 

223 evident, with the fraction of observations decreasing from more than 70% over the clear 

224 subtropical North Atlantic to about 35% over the cloudy equatorial Atlantic. After 

225 applying the data rejection described above, overall more than 50% of the observations 

226 remain. The most frequent rejections occur over the subtropical oceans, possibly due to 

227 the ambiguity between sea salt and dust aerosol over those regions 

228 It is noted that we applied Kaufman et a/.’s formula in a stricter way than what 

229 was originally done by Kaufman et al. [2005a and 2005b]. In Kaufman et al. [2005a] 
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230 when computing x smo ke, they set negative x sm oke values to zero, regardless of their 

231 magnitude. And in Kaufman et al. [2005b], when computing Xd us t, f is forced to be 

232 bounded within fd ust and f sm ok e by setting any f values outside the bounds to the limiting 

233 values, to guarantee Xdust is always positive. In our work, we include only those 

234 observations that fit the formula well. Nevertheless, it should be noted that although the 

235 data rejection we applied removes a large number of suspicious observations, doing so 

236 does not in itself guarantee the accuracy of the remaining applications, given that fd us t, 

237 fsmoke, and f mar are assumed constant, and x mar was estimated from an empirical formula. 

238 For the former factor, sensitivity tests will be performed to examine the extent to which 

239 our results are sensitive to the variations in FMF values. For the latter, the deviations of 

240 the actual x mar (mainly sea salt) compared to that empirically computed will be interpreted 

24 1 primarily as changes in dust amount, thus affecting the accuracy of Xdust- However, as will 

242 be presented next, MISR can distinguish dust from sea salt based on retrieved particle 

243 shape, thus Xd us t derived from MISR is not subject to the ambiguity between dust and sea 

244 salt that affects Kaufman’s method. In this sense the examination of MISR aerosol 

245 observations provides validation for the more extensive MODIS data set, in addition to 

246 offering actual results independent of MODIS. 

247 2.2. MISR 

248 The MISR instrument aboard the Terra satellite views Earth in four spectral bands 

249 in each of nine view angles spread out in the forward and aft directions along the flight 

250 path at 70.5°, 60.0°, 45.6°, 26.1°, and nadir. With a swath about 360 km wide, global 

251 coverage with MISR is acquired about once every 9 days at the equator and 2 days 

252 toward the poles. Data collected from all nine viewing angles provide independent 
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253 constraints on aerosol properties, separating mineral dust aerosols from other aerosol 

254 components operationally [Kahn et al., 2010]. 

255 Assuming the mineral dust is all non-spherical and the non-spherical part of AOT 

256 is all from dust, MISR Xd us t can be directly computed as the non-spherical fraction of the 

257 total AOT: 

258 Xdust — x x fnon-spherical (6). 

259 In MISR, smoke and marine aerosols (both sea salt and sulfate aerosols) 

260 contribute to the spherical part of the total AOT. (Refer to Table 3 in Kahn et al. [2001] 

261 for more details on shape categories of different aerosol components.) If x mar is again 

262 taken as computed from equation (5), x smo ke can be approximated as the difference 

263 between the total x and Xdust and x mar : 

264 X smo ke T — ^dust X ma r (7). 

265 However, it should be kept in mind that unlike MISR Xdust, the MISR x smo ke derived this 

266 way, and all the MODIS particle type distinctions, are not independent observations 

267 based on actual physical constraints, and their accuracy is subject to the large 

268 uncertainties due to the empirically calculated x mar . 

269 We use MISR-derived Level 3 daily x and non-spherical fraction at 0.558 pm on 

270 0.5° x 0.5° spatial grids during the same time period as the MODIS data. MISR x and Xd us t 

271 are interpolated onto 1° x 1° grids to compare with MODIS results. A number of 

272 validation studies have shown that overall, about 70% to 75% of MISR AOT retrievals 

273 fall within 0.05 or 20% X AOT, and about 50% to 55% are within 0.03 or 10% X AOT, 

274 except at sites where mixed dust and smoke are commonly found [Kahn et al., 2010], 
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275 Particle property validation suggests that expected MISR sensitivity to the spherical 

276 versus non-spherical particles is about 20% for AOT above 0. 15, and diminishes for mid- 

277 visible AOT below this value [Kahn et a/., 1997; Kalashnikova et ah, 2005; 

278 Kalashnikova and Kahn, 2006]. Thus we only use x > 0.15 to calculate idust from equation 

279 (6). With this x cutoff, it is found that more than 70% of the observations still remain for 

280 calculation over most of the tropical Atlantic (Figure 3b). Note generally MISR x 

281 retrievals are available for about 15% of the time, except over the convectively active 

282 regions (decreased to about 10%, Figure 3a). Note also that although it is never been 

283 explicitly addressed, MODIS sensitivity to particle properties, e.g. FMF, also diminishes 

284 at low AOT (implied in Figure 7 in Kahn et al. [2009]), and MISR actually has greater 

285 sensitivity at low AOT than MODIS due to the long atmospheric paths observed by its 

286 steeper-viewing cameras. 

287 Note that MISR also retrieves aerosol particle size infonnation, thus conceptually 

288 it is possible to follow the same method utilized for MODIS to separate the total AOT 

289 into specific types in MISR. However MISR categorizes the aerosol particles into three 

290 bins: “fine” (particle radii<0.35um), “medium” (radii between 0.35 and 0.7), and “large” 

291 (radii >0.7um) modes, instead of two bins as “fine” vs. “coarse” in MODIS. Thus 

292 applying Kaufman’s formula to MISR would require considerable additional work, but 

293 not necessarily lead to greater insight, because similar assumptions would be required to 

294 apply the size-discrimination method to MISR as to MODIS. Furthermore, the different 

295 radii range for MODIS and MISR “fine” mode would make it impossible to cross- 

296 validate the FMF values between these two datasets. Therefore, using MISR aerosol size 

297 information to separate different aerosol components following Kaufman’s method is 
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298 considered beyond the scope of current paper. 

299 3. Comparison of Aerosol Components between MODIS and MISR 

300 In this section, the AOTs of individual aerosol components over the tropical 

301 Atlantic Ocean are examined, and the MODIS and MISR results are compared. 

302 Climatological maps are examined first, in order to investigate how well the methods 

303 presented in Section 2 capture the basic features of the long-term mean. 

304 3.1. Aerosol Winter Climatology 

305 3.1.1. MODIS 

306 The MODIS climatological mean (2002-2009) boreal winter (November to April) 

307 aerosol maps over the tropical Atlantic Ocean for x, Xd us t, x smo k c , and x mar , as well as their 

308 percentages relative to x, are shown in Figure 4. The climatological mean x features a 

309 zonally oriented, optically thick aerosol plume centered at around 5°N-8°N stretching 

310 across the equatorial Atlantic Ocean. The magnitude and latitudinal extent are greatest 

311 over the eastern equatorial Atlantic along the west coast of Africa and gradually decrease 

312 westward toward the central and western equatorial Atlantic, as expected for an aerosol 

313 plume that originates in Africa (Figure 4a). Note that this map is similar to Figure lb 

314 (color shaded area) in Tian et al. [2011], but is not identical, due to the data rejection 

315 performed in this study. The magnitude of x in Figure 4a is about 80% of that in Tian et 

316 al. [2011] over the heavy aerosol loading region (x>0.15), and comparable over the clean 

317 subtropical Atlantic. Nevertheless, their overall patterns are similar, with spatial 

318 correlation as high as 0.98. The spatial pattern of climatological Xd us t resembles the x 

319 pattern greatly, with maximum Xd us t found along the west coast of Africa (about 0.35) and 

320 a gradual decrease toward the central and western equatorial Atlantic (Figure 4b). Over 
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321 the equatorial Atlantic Ocean, where aerosol loading is high (x>0.15), dust is the 

322 dominant aerosol component, contributing greater than 50%, and as much as 75%, to x 

323 (Figure 4c). 

324 Compared to Xd us t, Asmoke is significantly weaker. Plumes of fme-mode-dominant 

325 aerosol are found originating from the African (biomass burning smoke) as well as the 

326 South American continents (air pollution) (Figure 4d). A contribution of more than 20% 

327 to total x is found over the eastern tropical Atlantic, while less than 20% over the western 

328 part (Figure 4e). Figure 4f shows that x mar is very small (about 0.04) over the Atlantic 

329 intertropical convergence zone and the west coast of Africa because of the weak trade 

330 winds. Over the clean subtropical Atlantic, marine aerosol is the dominant component 

33 1 (>50%) due to lack of dust and smoke aerosols over these regions. 

332 These results indicate that the major aerosol plume over the equatorial Atlantic 

333 Ocean in Figure 4a is the dust originating in the Sahara desert and the Sahel region, with 

334 some contribution of biomass burning smoke originating in the Sahel and African 

335 savanna. These aerosol distributions are generally consistent with previous observational 

336 results [e.g., Husar et al., 1997; Kaufman et al., 2005b; Huang et al., 2010a], lending 

337 some confidence to the quantitative assessment of the relative dust, smoke and marine 

338 AOTs contributions to the total AOT from equations (3) and (4) from MODIS. 

339 3.1.2. MISR 

340 MISR winter climatologies of x, Xd us t, and x smo k e are shown in Figure 5. Recall that 

341 for these results only cases with x > 0.15 are used due to the fact that sensitivity to shape 

342 diminishes when aerosol concentration is low. Data rejection is also applied to MODIS x, 

343 too, as discussed in Section 2.1.3. Thus, we don’t expect Figure 5a to be the same as 
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344 Figure 4a. Nevertheless, we do find they have highly consistent spatial patterns and 

345 comparable magnitudes. The climatological MISR Xd us t map (Figure 5b) is also very 

346 similar to the MODIS idust (Figure 4b). The MISR x and Xdust are slightly larger than the 

347 MODIS counterparts partly due to MISR’s exclusion of low aerosol cases, but more 

348 important reasons will be addressed in Section 3.2. It is seen that MISR Xdust contributes 

349 more than 40%, up to more than 60%, to the total x for cases with x > 0.15, further 

350 confirming that dust is the dominant aerosol component over the equatorial Atlantic 

351 Ocean (Figure 5c). It is also noted that MISR has 10% to 15% lower dust fraction over 

352 the equatorial Atlantic. Again, this difference could be partly attributed to the x cutoff in 

353 MISR, but more importantly could result from both uncertainties involved in the 

354 derivation of MODIS and MISR Xdust, especially in MODIS, given the assumptions used. 

355 This point will be discussed again in Section 3.2. The climatological MISR x smo k e , 

356 calculated as described in Section 2.2, has a 10% larger contribution to total x than 

357 MODIS x sm oke over the equatorial Atlantic (Figure 5d and 5e compared to Figure 4d and 

358 4e). 

359 The overall consistency between the MODIS and MISR aerosol climatology gives 

360 us confidence in applying two independent satellite datasets and methods to derive the 

361 dust and smoke aerosols over the tropical Atlantic Ocean. However, the consistency in 

362 climatology does not guarantee their consistency on shorter time scales, for example, on 

363 the intra-seasonal time scale of importance here. Therefore, in next subsection we directly 

364 compare the coincident daily MODIS and MISR Xdust over the tropical Atlantic Ocean. 

365 3.2. Comparison of Coincident MODIS and MISR Dust Aerosols 

366 In order to compare the daily Xd us t derived from MODIS and MISR, the 
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367 correlation between them for all their coincident days during 2002-2009 over the tropical 

368 Atlantic is calculated. Here, the term “coincident” simply means the MODIS and MISR 

369 observations fall within a same grid box and on a same day, which is less strict than its 

370 usual implication used by the satellite community. We use the whole years for 2002-2009 

371 instead of winters only in order to obtain as many coincident days as possible. The 

372 correlation is only calculated when there are at least five coincident days at a grid point. 

373 Figure 6a shows that overall the MODIS and MISR dust AOTs are well 

374 correlated. The correlation is systematically higher in the north Atlantic region than that 

375 in the south, with more coincident days over the north Atlantic (Figure 6b). Over the 

376 North Atlantic, the correlation coefficient is typically around 0.7 or larger, whereas over 

377 the South Atlantic, the correlation is typically less than 0.5 with quite a few spots less 

378 than 0.2, or even negative. Generally, the correlation is higher when it is closer to the 

379 African continent, and decreases gradually as the dust is transported away from the 

380 source of the dust load. These results suggest that the MODIS and MISR derived dust 

381 aerosols agree with each other quite well over heavy dust load regions, whereas they are 

382 less consistent over the regions with less frequent dust occurrence or small dust aerosol 

383 concentration. 

384 Figure 6a shows that overall the MODIS and MISR Xd us t are highly correlated, 

385 however the correlation coefficient does not provide information on the Xd us t magnitude. 

386 Therefore, the time series of MODIS and MISR Xd us t as the function of their coincident 

387 days at a representative grid point (29.5°W, 19.5°N), as well as averaged within the 5° x 

388 5° and 10° x 10° gird boxes centered at this point, are shown (Figure 6c-e). Again it is 

389 found that the MODIS and MISR Xd us t are highly correlated on the daily basis, and it also 
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390 reveals that Xdust is systematically greater for MISR than MODIS. 

391 MISR idust is larger than MODIS idust over most of the tropical Atlantic when 

392 averaged for coincident days (about 0.03 larger averaged over the basin), except in some 

393 regions north of the Equator (Figure 7b). This difference can be traced back to the 

394 difference between MISR and MODIS total AOT (Figure 7a). Over almost the entire 

395 tropical Atlantic, coincident MISR AOTs are systematically larger than the MODIS ones 

396 (up to 0.08, and about 0.04 when averaged for the Atlantic basin). This is consistent with 

397 previous studies in which MISR AOT is found to be generally larger than MODIS AOT 

398 over water [e.g., Abdou et al., 2005; Kahn et al., 2010]. Further examination of 

399 coincident MODIS and MISR x and Xdust binned against the MISR x reveals that larger 

400 MISR x is found for the entire x spectrum, whereas MISR Xd us t is larger (smaller) when 

401 the aerosol concentration is relatively low (high) (Figure 7c). Note that lower aerosol 

402 concentrations are overwhelmingly more frequent (black line in Figure 7c), thus the 

403 averaged MISR Xd us t is larger than the averaged MODIS Xd us t- As seen above, the Xd us t 

404 difference between MODIS and MISR can be traced back to the x difference between 

405 them. However, this is not the only reason, as, unlike MISR x, MISR Xd us t is not 

406 systematically larger than MODIS Xd us t- The uncertainties involved in the derivation of 

407 both MODIS and MISR Xd us t inevitably lead to their differences too. However it should be 

408 noted that although both methods have limitations, MISR discrimination of aerosol type 

409 is much more robust than MODIS discrimination. The derivation of MISR Xdust is based 

410 on actual retrieved particle property information, whereas for MODIS, the accuracy of 

411 the formula is subject to the use of constant FMF values and empirically calculated x mar . 

412 Despite more physically robust separation of Xdust from x in MISR, it has much less 
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413 frequent sampling compared to MODIS, thus it is necessary to examine both datasets. 

414 Results based on these two complementary satellite datasets will provide more solid 

415 insight to the characteristics of the MJO-related dust and to some extent smoke 

416 variabilities. 

417 4. MJO-related Atlantic Dust and Smoke AOT Anomalies 

418 4.1. MJO Analysis Methodology 

419 For the MJO analysis and composite procedure, we use the multivariate empirical 

420 orthogonal function (EOF) method introduced by Wheeler and Hendon [2004] and 

421 adopted widely by the MJO community [e.g., Waliser et al., 2009]. Briefly, the intra- 

422 seasonal anomalies of daily AOT are obtained by removing the climatological-mean 

423 seasonal cycle and filtering via a 30-90-day band pass filter. Then, a composite MJO 

424 cycle (8 phases) is calculated by averaging the daily anomalies that occur within each 

425 phase of the MJO cycle. The MJO phase for each day is determined by the Real-time 

426 Multivariate MJO (RMM) index (a pair of PC time series called RMM1 and RMM2). 

427 Only days having strong MJO activity (RMM1 2 +RMM2 2 >=1) are considered. 

428 Figure 8 shows the number of strong MJO days, or events 

429 (RMM1 2 +RMM2 2 >=1), used for the 8-phase MJO cycle composite for both MODIS and 

430 MISR during 2002-2009 boreal winters. The number of total strong MJO events for each 

431 MJO phase during this period is also shown at the right upper comer of each panel. For 

432 both instruments, the number of strong MJO events used for the composite is much less 

433 than the actual number of total events because of satellite retrieval sampling issues and 

434 the data rejection applied. Generally, the number of MODIS events ranges from about 10 

435 to 45, about as three times more than MISR, due to the much wider MODIS swath. 
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436 4.2. MODIS and MISR 

437 The 8-phase MJO composite maps of MODIS total x, idust and x smo ke anomalies are 

438 shown in Figure 9a, 10a and 1 la, respectively. The MISR counterparts will be discussed 

439 and compared with MODIS results later in Section 4.4. Comparing Figure 9a and 10a, it 

440 is seen that MODIS x and Xdust have very similar temporal evolution and spatial patterns, 

44 1 except that Xdust anomalies are slightly smaller. The spatial correlation between x and Xdust 

442 anomalies for 8 MJO phases is 0.89. For both x and Xdust, strong negative anomalies (as 

443 large as about -0.04) are found over the entire equatorial Atlantic for MJO phases 1, 2, 3, 

444 and 8. In contrast, strong positive anomalies (up to 0.04) are found over the equatorial 

445 Atlantic for MJO phases 5-6. For MJO phase 4, strong positive anomalies occur to the 

446 north of the equator, whereas negative x anomalies are found to the south. The converse 

447 is true for the MJO phase 7. The MJO composite maps of total x greatly resemble those 

448 shown by Tian et a/. [2011], with spatial correlation 0.52 (significant at the 99.9% level), 

449 and have very similar magnitudes. This resemblance suggests that the MJO-related total 

450 AOT anomaly patterns are robust and not sensitive to the data sampling (fewer samples 

45 1 are used in this study). The MJO-related x smo k e anomalies in MODIS are very weak, rarely 

452 exceeding 0.01(Figure 11a). The spatial correlation between x smo k e and x is only about 

453 0.21. The MJO composite maps of x mar anomalies are not shown since x mar is linearly 

454 dependent on the surface wind speed, thus the MJO-related x mar anomaly pattern in fact 

455 reflect the wind anomalies associated with the MJO, which has been examined in Tian et 

456 al. [2011], Furthermore, it is found that the magnitudes of the x mar anomalies are 

457 negligible: the strongest negative/positive anomalies are about -0.004/0.004 (figure not 

458 shown). 
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459 The 8-phase MJO composite maps of MISR total x, idust and i smo ke anomalies are 

460 shown in Figures 9b, 10b and lib. Again, it is found that the Xdust anomalies are 

461 significantly larger than the x smo k e anomalies. The Xdust anomalies have very similar 

462 patterns to those of total x, with slightly smaller magnitude, whereas the x smo ke anomalies 

463 are very small and noisy. Further comparison between the MISR and MODIS x and Xdust 

464 results (Figures 9 and 10) indicates that overall they exhibit very similar temporal 

465 evolution as well as anomaly patterns despite the systematically larger MISR anomalies 

466 compared to those of MODIS, which is likely the outcome of the systematically larger 

467 MISR x and Xdust retrievals compared to MODIS over the tropical Atlantic Ocean, as 

468 discussed earlier. The rejection of all x < 0. 15 data in MISR could contribute to the above 

469 difference, too; however, MISR x anomalies with no data rejection don’t show evident 

470 differences (Figure not shown). Besides the magnitude difference, MISR x and Xdust 

471 anomalies are also noisier than the MODIS anomalies due to lower sampling. 

472 Nevertheless, the consistency between the MODIS and MISR results shown in Figures 9- 

473 11 demonstrate that dust is the dominant aerosol component on the intra-seasonal time 

474 scale, and the MJO-related dust anomalies are robust, as seen from two independent sets 

475 of satellite observations with different strengths and limitations. 

476 4.3. Sensitivity of MODIS Results 

477 Kaufman et al.’s method to compute Xdust and x smo k e using MODIS x and f based 

478 on equations (3) and (4) is straightforward; however, large uncertainties in the computed 

479 Xdust and x smo k e are expected for at least several reasons: the uncertainties in the MODIS x 

480 and f retrievals, the uncertainties in empirically computed x mar , and most critically, the 

48 1 uncertainties resulted from assuming constant f mar , f dust , and f smo k e . In this subsection, we 
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482 will examine the sensitivity of our results to the FMF values. We first perturb the FMF 

483 values used in this study (see Section 2.1.1) by increasing/decreasing one of them at a 

484 time by its uncertainty range while keeping the other two unchanged. These sensitivity 

485 test cases are denoted group 1. We then test using FMF values derived by three other 

486 studies [Jones and Christopher, 2007; 2011; Yu et ah, 2009], which are denoted group 2. 

487 These studies have attempted to calibrate Kaufman's technique, and re-derived the FMFs 

488 using either updated MODIS datasets, or aerosol observations from other satellite 

489 datasets, or using the GOCART (Goddard Chemistry Aerosol Radiation and Transport) 

490 model to locate regions where a single aerosol component dominates. The FMF values 

491 they obtained are generally consistent with what is used in this study: f mar , fdust, and f smo ke 

492 are 0.25, 0.44, and 0.83, respectively in the study by Jones and Christopher [2007], and 

493 0.31, 0.49, and 0.78 in their 2011 study, and 0.45, 0.37, and 0.90 in Yu et al. [2009]. 

494 FMFs do vary, and these studies find that they vary considerably depending on region 

495 and season. 

496 Figure 12a shows the composite 8-phase MJO cycles of Xd us t anomalies averaged 

497 over the equatorial Atlantic (30°W-15°W, EQ-15°N) based on different sensitivity test 

498 cases. The 30°W-15°W, EQ-15°N box is chosen to represent the region with strong 

499 intra-seasonal aerosol variations. Overall, the colored lines (9 sensitivity tests) cluster 

500 around the solid black curve (control case), and the MJO cycles of Xd us t anomalies based 

501 on different sensitivity tests show a coherent evolution. This suggests that the MJO 

502 associated Xdust anomalies over the tropical Atlantic are quite robust despite the 

503 uncertainties in using constant FMFs. Nevertheless, the lines spread in some phases. The 

504 result based on Yu et al. [2009] is most different from other cases, probably because in 
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505 the other studies, the FMF values follow the sequence f mar < fdust < fsmoke despite the 

506 deviations, whereas the order of f mar and fdust is reversed in Yu et a/.’s study [2009]. The 

507 spread in group 2 is naturally larger than that in group 1. Furthennore, the idust anomalies 

508 are more sensitive to fd USt , as indicated by the larger deviations of the dashed cyan and 

509 blue curves to the control case in phase 2, 3 and 6. 

510 Similarly, Figure 12b shows the sensitivity test results for T smo k e anomalies. The 

5 1 1 spread of x sm ok e anomalies is quite large compared to their magnitude. Nevertheless, the 

512 overall MJO cycle of x smo ke anomalies is consistent among the different cases and their 

513 magnitudes are much smaller than that of the Xdust anomalies. 

5 14 5. Summary and Conclusions 

515 Previous studies [Tian et al., 2008; 2011] found significant intra-seasonal 

516 variability related to the MJO in the total column AOT over the tropical Atlantic region. 

517 Aerosol over the tropical Atlantic is primarily a mixture of mineral dust, biomass burning 

518 smoke, and marine aerosol. Given the different roles in the radiative forcing and cloud 

519 formation processes played by these three aerosol types, as well as the potential 

520 predictability of the MJO extending to 2-4 weeks, it is of great interest to further examine 

521 the MJO-related variability for individual aerosol types, especially dust and smoke 

522 aerosols. 

523 Daily MODIS/Aqua total AOT and FMF measurements are used to derive daily 

524 Xdust and x smo k e following Kaufman et a/.’s method [2005a, b]. This method contains 

525 considerable uncertainties, due to the assumption of constant f dust , f smo ke , fmar as well as 

526 the empirical calculation of x mar . Strict data rejection has been applied in order to use 

527 Kaufman’s formula in a safer way. With MISR’s sensitivity to aerosol particle shape, 


23 



528 dust and smoke aerosols can be distinguished using the MISR aerosol non-spherical 

529 fraction. Because MISR sensitivity to shape (and MODIS sensitivity to FMF) diminishes 

530 when aerosol concentration is low, only x >= 0. 15 data are used to compute idust and Tsmoke 

531 for MISR. The examination of both datasets is necessary, and results from MODIS and 

532 MISR are complementary: MODIS provides better spatial coverage, and MISR dust is 

533 derived from actual aerosol property retrieval rather than assumed aerosol-type-specific 

534 FMF. 

535 Both MODIS and MISR show a very similar dust and smoke winter climatology. 

536 Xdust is found to be the dominant aerosol component over the tropical Atlantic. It is largest 

537 over the eastern equatorial Atlantic (about 0.35) and gradually decreases toward the west. 

538 x sm oke is significantly smaller than Xdust, Over the equatorial Atlantic with strong total 

539 aerosol loading, the contribution of Xdust to total x ranges from 40% to 70%, in contrast to 

540 about 25% is attributed to x smo k e and less than 20% from x mar - The daily MODIS and 

541 MISR Xdust are overall highly correlated, with the correlation coefficients typically about 

542 0.7 over the North Atlantic, but much smaller or even negative over the South Atlantic. 

543 MISR Xdust is found to be systematically greater than the coincident MODIS Xdust, and this 

544 difference can be traced to the AOT difference between them. The consistency of the 

545 MODIS and MISR dust and smoke aerosol climatology and daily variations give us 

546 confidence to use these two data sets to investigate the relative contribution of dust and 

547 smoke aerosols to the total AOT variation associated with the MJO. However, the 

548 identification of smoke is much less certain than that of dust, because discrimination 

549 among fine -mode sea salt, sulfate, and smoke particles depends on assumptions for both 

550 MODIS and MISR, whereas the MISR dust discrimination is based on retrieved particle 
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551 


shape. 


552 For MODIS, the MJO composite maps of idust anomalies are very similar to those 

553 of x anomalies, and are of comparable magnitude. Furthermore, the variance of idust 

554 anomalies on the full intra-seasonal time scale is found to be comparable or even bigger 

555 than that of the x anomalies. In contrast, the MJO-related x smo ke anomalies are rather 

556 small, barely exceeding 0.01, and the x mar anomalies are negligible. The sensitivity study 

557 further shows that the MJO-related Xdust and x smo ke anomalies are quite robust, even when 

558 we perturb the FMFs by their uncertainty ranges or use different sets of FMFs from 

559 several independent studies. 

560 Similarly, MISR also shows that the MJO composite maps of Xdust anomalies are 

561 very similar to those of x anomalies, while the MJO-related x smo ke anomalies are rather 

562 small. The composite MJO cycle of Xdust anomalies from MISR over the tropical Atlantic 

563 Ocean is consistent with the MODIS results although the anomalies are much noisier due 

564 to its less frequent sampling. The magnitude of MISR anomalies is again found to be 

565 systematically larger than that of MODIS. The MJO-related x smo k e anomalies in MISR are 

566 overall slightly larger than for MODIS, but still much smaller compared to MISR Xd us t 

567 anomalies. 

568 The consistency between the MODIS and MISR Xdust and tsmoke anomalies in 

569 terms of the evolution of the MJO cycle as well as the spatial pattern of anomalies 

570 suggests that dust aerosol is the dominant component on the intra-seasonal time scale 

57 1 over the tropical Atlantic Ocean. 

572 The observational results obtained from two complementary satellite datasets can 

573 be used to evaluate chemical transport models and help in model development. 
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574 Furthermore, given the potential predictability of the MJO, this study implies the MJO 

575 modulation of Atlantic dust aerosol concentration may be predictable with lead times of 

576 2-4 weeks, which in turn may lend guidance to many other dust-related phenomena, such 

577 as the air quality, dust storm activity, and ocean nutrient deposition over the Atlantic 

578 Ocean. 
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742 Figure Captions 

743 Figure 1: Calculated a) Xd us t and b) x smo ke as a function of total x and f according to 

744 equations (3) and (4). Here, f mar =0.3, fdust=0.5, f sm oke = 0.9, and x mar =0.06, as indicated by 

745 four white straight lines. The white dashed lines indicate Xd us t =_ 0.03 and x smo ke =_ 0.03. 

746 Contour interval: 0.1. 

747 Figure 2: a) The count of paired MODIS x and f measurements with respect to x and f 

748 over the tropical Atlantic (20°S-30°N, 60°W-20°E) for 2002-2009 boreal winters. 

749 Superimposed solid lines are the lines of Xd U st=~0.03 and x sm oke =_ 0.03 in Figure 1. b) The 

750 percentage of days with MODIS x and f measurements during 2002-2009 boreal winters 

75 1 (1269 days in total), c) The percentage of x and f measurements used to calculate Xd us t and 

752 x sm oke relative to the total number of measurements. 

753 Figure 3: a) The percentage of days with MISR x measurements during 2002-2009 

754 boreal winters (1269 days in total), b) The percentage of MISR x > 0.15 relative to total 

755 number of x measurements. 

756 Figure 4: MODIS climatological mean (2002-2009) boreal winter a) total x, b) Xd us t, d) 

757 x smo ke, and f) x ma r over the tropical Atlantic Ocean, as well as their percentages to x (c- 

758 Xdust, o -x S moke, and g-Xmar)- 9-point spatial smoothing is applied. 

759 Figure 5: Same as Figure 3a-3e but for MISR, and only cases with MISR x > 0.15 used. 

760 Figure 6: a) Correlation coefficients between coincident MODIS Xdust and MISR Xd us t 

761 during 2002-2009 over the tropical Atlantic and b) the number of coincident days. Time 

762 series of MODIS (blue) and MISR (red) Xd us t as the function of the coincident day c) at 

763 the point of 29.5°W, 19.5°N, averaged on the d) 5° x 5° and e) 10° x 10° boxes centered 

764 at 29.5°W, 19.5°N. 
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765 Figure 7: a) MISR and MODIS x difference averaged over coincident days during 2002- 

766 2009. b) Same as a), but for idust- c) Number of coincident MODIS and MISR 

767 observations (black line), averaged MODIS x (thick blue line), MISR x (thick red line), 

768 MODIS Xdust (thin blue line), and MISR Xdust (thin red line) over 20°S-30°N, 60°W-20°E 

769 during 2002-2009 as a function of binned MISR x (binned by every 0.01). 

770 Figure 8: The number of strong MJO events used for MJO composite for each phase of 

771 the MJO cycle for a) MODIS and b) MISR during the 2002-2009 boreal winters. The 

772 total number of strong MJO events during this period is indicated at the right upper 

773 comer of each panel. 

774 Figure 9: MJO composite maps of total x anomalies (multiplied by 100) for a) MODIS 

775 and b) MISR over the tropical Atlantic Ocean. 9-point spatial smoothing is applied. 

776 Figure 10: Same as Figure 9, but for Xdust- 

777 Figure 11: Same as Figure 9, but for x sm oke- 

778 Figure 12: The composite MJO cycle of MODIS a) Xd us t and b) x smo k e anomalies averaged 

779 over 30°W-15°W, EQ-15°N for the control case (black) and 9 sensitivity test cases 

780 (color) with f mar , fd ust , and f sm oke of each case indicated in the legend. 

781 
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Figure 1: Calculated a) Td us t and b) T sm oke as a function of total t and f according to 
equations (3) and (4). Here, f mar =0.3, fd us t=0.5, f sm oke = 0.9, and x mar = 0.06, as indicated by 
four white straight lines. The white dashed lines indicate Tdust =_ 0.03 and T sm oke =_ 0.03. 
Contour interval: 0.1. 
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789 Figure 2: a) The count of paired MODIS x and f measurements with respect to x and f 

790 over the tropical Atlantic (20°S-30°N, 60°W-20°E) for 2002-2009 boreal winters. 

791 Superimposed solid lines are the lines of Xdust=~0.03 and x sm oke =- 0.03 in Figure 1. b) The 

792 percentage of days with MODIS x and f measurements during 2002-2009 boreal winters 

793 (1269 days in total), c) The percentage of x and f measurements used to calculate Xd us t and 

794 x smo ke relative to the total number of measurements. 
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Figure 3: a) The percentage of days with MISR x measurements during 2002-2009 
boreal winters (1269 days in total), b) The percentage of MISR x > 0. 15 relative to total 
number of x measurements. 
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802 

803 Figure 4: MODIS climatological mean (2002-2009) boreal winter a) total x, b) Xd us t, d) 

804 x smo ke? and f) x mar over the tropical Atlantic Ocean, as well as their percentages to x (c- 

805 Xdust, e-Xsmoke, and g-x mar ). 9-point spatial smoothing is applied. 


39 


806 



40 


809 





810 Figure 6: a) Correlation coefficients between coincident MODIS Xd us t and MISR Xd us t 

811 during 2002-2009 over the tropical Atlantic and b) the number of coincident days. Time 

812 series of MODIS (blue) and MISR (red) Xdust as the function of the coincident day c) at 

813 the point of 29.5°W, 19.5°N, averaged on the d) 5° x 5° and e) 10° x 10° boxes centered 

814 at 29.5°W, 19.5°N. 
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816 Figure 7: a) MISR and MODIS t difference averaged over coincident days during 2002- 

817 2009. b) Same as a), but for Xdust- c) Number of coincident MODIS and MISR 

818 observations (black line), averaged MODIS t (thick blue line), MISR x (thick red line), 

819 MODIS idust (thin blue line), and MISR idust (thin red line) over 20°S-30°N, 60°W-20°E 

820 during 2002-2009 as a function of binned MISR x (binned by every 0.01). 
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822 Figure 8: The number of strong MJO events used for MJO composite for each phase of 

823 the MJO cycle for a) MODIS and b) MISR during the 2002-2009 boreal winters. The 

824 total number of strong MJO events during this period is indicated at the right upper 

825 comer of each panel. 
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827 Figure 9: MJO composite maps of total t anomalies (multiplied by 100) for a) MODIS 

828 and b) MISR over the tropical Atlantic Ocean. 9-point spatial smoothing is applied. 
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83 1 Figure 10: Same as Figure 9, but for Xd us t- 
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835 Figure 11: Same as Figure 9, but for T smo k e - 
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Figure 12: The composite MJO cycle of MODIS a) idust and b) T smo ke anomalies averaged 
over 30°W-15°W, EQ-15°N for the control case (black) and 9 sensitivity test cases 
(color) with f mar , fdust, and f sm oke of each case indicated in the legend. 
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